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ABSTRACT 



Aims. We study the morphology of a set of superclusters drawn from the SDSS DRV. 

Methods. We calculate the luminosity density field to determine superclusters from a flux-limited sample of galaxies from SDSS 
DR7, and select superclusters with 300 and more galaxies for our study. We characterise the morphology of superclusters using the 
fourth Minkowski functional V3, the morphological signature (the curve in the shapefinder's Ki-K^ plane) and the shape parameter 
(the ratio of the shapefinders K1/K2). We investigate the supercluster sample using multidimensional normal mixture modelling. We 
use Abell clusters to identify our superclusters with known superclusters and to study the large-scale distribution of superclusters. 
Results. The superclusters in our sample form three chains of superclusters; one of them is the Sloan Great Wall. Most superclusters 
have filament-like overall shapes. Superclusters can be divided into two sets; more elongated superclusters are more luminous, richer, 
have larger diameters, and a more complex fine structure than less elongated superclusters. The fine structure of superclusters can be 
divided into four main morphological types: spiders, multispiders, filaments, and multibranching filaments. We present the 2D and 
3D distribution of galaxies and rich groups, the fourth Minkowski functional, and the morphological signature for all superclusters. 
Conclusions. Widely different morphologies of superclusters show that their evolution has been dissimilar. A study of a larger sample 
of superclusters from observations and simulations is needed to understand the morphological variety of superclusters and the possible 
connection between the morphology of superclusters and their large-scale environment. 
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1. Introduction 

The most remarkable feature of the megaparsec-scale mat- 
ter distribution in the Universe is the presence of the cosmic 
web - the network of galaxies, groups, and clusters , connected 
by filaments ( Joeveeretal. 1978; Gregory & Thom pson! 1 19781 
IZeldovich et a l. 1982; de Lapparentet al. 1986). The formation 
of a web of galaxies and systems of galaxies is predicted in any 
physically motivate d theory of the for mation of structure in the 
Universe (see, e.g.. iBond et al.lll996h . In this scenario galax- 
ies and galaxy systems form because of initial density pertur- 
bations on different scales. Perturbations on a scale of about 
100 /i"' Mpc (Ho = 100/zkms^'Mpc"') give rise to the largest 
systems of galaxies - rich superclusters. At larger scales dynam- 
ical evolution proceeds at a slower rate and superclusters have 
retained the memory of the initial conditions of their forma- 
tion and of the early evolution of structure (lEinasto et al.lll980l: 
IZeldovich et alJll982tlKoftnan et alJll987h . 

Numerical simulations show that high-density peaks in the 
density distribution (the seeds of supercluster cores) are seen 
already at very ea rly stages of the formation and evolution of 
structure (lEinastoll2010l) . These are the locations of the forma- 
tion of_tfie_^rs^_obiectsjiithe^^ (e.g. IVenemans et aP 
I2004t iMobasher et al.1 1200I lOuchi et al.1 [20051) 7 Observations 
have already found superclusters at high redshifts (Nakata et al. 
[2005; Swinbank et al. 2007: G al et al. 2008: Tanaka et aL2009) . 



Recently, the XMM-Newton satellite follow-up for the vali- 
dation of Planck cluster candidates led to the discovery of 
two massive, previously unknow n superclusters of galaxies 
jPlanck Collaboration et al.l|201 1|) . These are likely the first su- 
perclusters discovered through the Sunyaev-Zeldovich effect. 

Superclusters are important tracers of dark and bar^- 
onic matter in the Universe (Zappacosta et al. 20051 
Genova-Santos et al. 2005; Hevmans et al. 2008; Buote et al] 
2009; Pad illa-Torres et al. 2009; Schirm er et a l. 2011). Studies 
of superclusters and of the supercluster-void network have 
demonstrated the presence of a chara cteristic scale in th e 
distribution of rich superclusters (Einasto ~et al.l Il994[ Il997ah . 
This was probably an early hint of baryon acoustic oscillations 
(,Hiitsiii201Q) . 

To search for superclusters and to understand their prop- 
erties, we need to know how to identify them and how to 
quantify their properties (.Bond et al.l l2010i) . Several m ethods 



have b een proposed to study the cosmic web ( Bharadw ai et al. 



a et aljE 
bie et al] 



20101: lAragon-Calvo et al.1 12010c iSousbie 
120111 and references therein). One ap 
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200l; -Stoic 
2011; Sousbie 

proach is to determine cosmic structures (in our study - 
superclusters of galaxies) using the density field and to 
study their morphology with Minkowski functionals and 
shapefinders (Schmalzing & Buchert 1997; Sathvaprakash et ajJ 
1998; Basilakos 2003; Sheth et al. 2003; Shandarin et al. 200^ 
rEinastoetal]|2007dl and references therein). fOortI (il983l) gave 
a review of the early studies of superclusters. Supercluster 
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catalogues compile d based on data ab o ut clusters of g alaxie s 
were published by Zucca et al.' (1993); Einasto et alJ (Il994l) : 
[Kalinkov & Kuneva (1995); Einasto et al. (2001). Recent deep 
surve y s of galaxies (as 2dF GRS and SDSS, see ICoUess et al.l 
120031; lAbazaiian et al.l 12009.) introduced a new era in the 
studies of the large-scale structure of the Universe, where 
systems of galaxies can be studied in unprecedended de- 
tail. A number of supercluster catalogues have bee n com- 
piled with these da ta, (Basilakos 2003; Einasto e t al. 2003at 



Erdogdu et al."2004!: 'Einasto et al..,2006> .2007^ Liivamagi et all 
2010 ; Luparello et al. 2011, and references therein). 

The overall morphology of superclusters has been studied 
by several authors (Kolokotronis et al. 2002; Basi lakos 2 003; 



ICosta-Duarte et al.ll201 ih . we refer to lEinasto eTaLr(l2007al) for 
a review of the properties of superclusters. The shapes and 
sizes of superclusters can be used to compare the observed su- 
perclusters with tho se obtained from cosmologic al simulations 
Kolokotronis et alj [2002; Einasto et al.1 l2007ah . fNichol et al 



2006h showed that the higher order correlation functions of 



the 2dFGRS do not agree with those found in numerical sim- 
ulations (but this fact can be explained by non-Gaussian ini- 
tial density fields. lOaztanaga & Maehoene n 1996). This discrep- 
ancy may be caused by the unusual morphology of one of the 
richest superclusters in the 2dFGRS, the supercluster SCI 126 
dEinasto et al."2007d', '2008") from the catalogue of superclusters 
bylEinasto et al. (2001, hereafter EOl). The morphology of su- 
perclusters may b e used to distinguish be tween different cosmo- 
logical models (iKolokotronis et afl2002h . 

Superclusters contain structures with a wide range of 
densities, from high-density cores of rich clusters to low- 
density filaments between clusters and groups. This makes 
them ideal laboratories to study processes that affect the evo- 
lution of galaxies, groups, and clusters of galaxies. A num- 
ber of studies have already shown that the supercluster envi- 
ronment affects t he properties of galaxies, groups, and clus- 
ters located there (lEinasto et al.ll2003k ■PUonis.,2004i: IWolf et al 



2005 1 ; [Haines et al J l2006t lEinasto et all l2007d; iPorteretal 
20081; iTempel et al.' '2009^, 'Fleenor & J ohnston-HollittI 1201 Ot 
Tempel et al. 2011; Einasto et al. 2011). Other evidence about 
the influence of the large-scale environment of galaxies on their 
properties co mes from the study of the properties of galaxies 
in void walls dCeccarelli et al.ll2008l) and from the study of the 
large-scale environment of quasars (' Lietzen et al.ll2009l) . A de- 
tailed information on the morphology of superclusters is needed 
to find out whether that may also be an important environmental 
factor in shaping the propert ies of galaxies and groups of galax- 
ies in superclusters (see a lso lAragon-Calvo et al. mm . 

lEinasto et alj (l2007dl) showed that the morphology of a typ- 
ical poor supercluster can be described as a "spider" - a sys- 
tem of several filaments growing from one concentration cen- 
tre (a rich cluster). The Local Supercluster is an example of a 
typical "spider'. Rich superclusters can be described as "mul- 
tispiders", where several high-density clumps are connected by 
lower-density filaments. One very rich supercluster (SCI 126) 
was described as a multibranching filament that consists of a 
rich filament of a quite uniform high density with poorer fila- 
ments for branches. This work concerned only a few of the rich- 
est supercluster while the present analysis enables us to test the 
classification of the morphology of superclusters. 

The goal of the present paper is to study in detail the 
morphology of a large sample of superclusters drawn from 
the 7th data release of the Sloan Digital Sky Survey (SDSS). 
The superclusters are determined using the global luminosity 
field and their morphologies are quantified with the Minkowski 



functionals and shapefinders. The superclusters are divided 
into two sets using multidimensional normal mixture mod- 
elling, applying the Mclust package for clustering and clas- 
sification dFralev & Raftervl 120061) from R, an open-source 
free statistical environmen t developed under the GNU GPL 
dlhaka & Gentleman! 1 19961 http://www.r- project.org). 
The cluster membership of superclusters and their large-scale 
distribution is analysed, and a short description of individ- 
ual superclusters is given. The data are described in Sect. |2] 
the methods in Sect. |3l and the results in Sect. |4| In Sect. |5] 
the selection effects in our sample are discussed, and a 
comparison with other studies is given. The study is sum- 
marised in Sect. |6] Interactive 3D models of the richest su- 
perclusters in our sample can be found on our web page: 
http : //www . aai . ee/~maret/SDSSsclmorph . html. 

We assume the standard cosmological parameters: the 
Hubble parameter Hq - 100 /z km s ' Mpc the matter den- 
sity Qtn = 0.27, and the dark energy density Qa = 0.73. 



2. Data 

We selected the MAIN galaxy sample of the 7th data release of 
the Sloan Digital Sky Survey dAdelman-McCarthv et aP 120081: 



lAbazaiian et al.ll2009l) with the apparent r magnitudes 12.5 
r < \1 .11, excluding duplicat e entries. The sample is described 
in detail in iTago et al.l d2010l) . hereafter TIO. We corrected the 
redshifts of galaxies for the motion relative to the CMB and 
computed the co-moving distances dMartmez&Saad 120021) of 
galaxies. 

The absolute magnitudes of galaxies are determined in 
the r-band M, with A:-correction for the SDSS galaxies cal- 
culated with the KCORRECT algorithm (< Blanton et a"l]|2003al ; 
Blanton & Roweis 2007). In addition, we applied the evo- 
lution correct i ons, using the luminosity evolution model of 
iBlanton et alj d2003b.) . The magnitudes correspond to the rest- 
frame at the redshift z = 0. 

The first step is to determine groups and clusters of galaxies 
with the friends- of-friends algorithm, where a galaxy belongs to 
a group of galaxies if this galaxy has at least one group member 
galaxy closer than the selected linking length. The linking length 
along with the distance was increased, to take into account se- 
lection effects, when constructing a group catalogue for a flux- 
limited sample. As a result, the maximum sizes and velocity dis- 
persions of groups are similar at all distances. For details and for 
the group catalogue we refer the reader to TIO. Q 

To determine the luminosities of groups and to calculate the 
luminosity density field we have also to correct for the luminosi- 
ties of galaxies that lie outside of the survey magnitude range. 
The calculation of luminosities is describe d in Appendix[Al (de- 
tails of this calculation are given also in Tempe l et al.ll20I ll) . 

In the final flux-limited group catalogue the richness of 
groups decreases rapidly at distances D > 320 Mpc because 
of selection effects. This effect is seen in Fig. [1] At small dis- 
tances, D < 70 /z"' Mpc luminosity weights are large owing to 
the absence of very bright galaxies. Therefore we chose for the 
present analysisa subsample of galaxies and galaxy systems in 
the distance interval 90 /; ' Mpc < D < 320 h'^ Mpc where the 
selection effects are small. 



' The TIO group catalogue is available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
or via http ://cdsweb .u-strasbg . fr/cgi-bin/qcat?J/A+A/ 
514/A1I82. 
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Table 1. Data of the superclusters. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


5 CliD 


IDeoi 


ID 


^gal 


^30 


"^peak 
Mpc//-! 


hot 


Diam. 
Mpc//i 


'^magV 
mag 


^galV 


Vs.max 




Ki 


KJK2 


1 


162 


239+027+009 


1038 


4 


264.5 


1591.5 


50.3 


-19.70 


916 


2 


0.077 


0.163 


0.48 


10 


160 


239+016+003 


1463 


2 


111.9 


680.2 


22.8 


-17.50 


1455 


1 


0.026 


0.024 


_ 


11 


154 


227+006+007 


1222 


5 


235.0 


1476.0 


35.4 


-19.30 


1180 


3 


0.051 


0.056 


0.92 


24 


111 


184+003+007 


1469 


4 


230.4 


1768.2 


56.4 


-19.25 


1419 


6 


0.077 


0.163 


0.48 


55 


111 


173+014+008 


1306 


4 


242.0 


1773.0 


50.3 


-19.55 


1155 


5 


0.084 


0.198 


0.43 


60 


160 


247+040+002 


1335 


3 


92.0 


527.4 


21.2 


-17.30 


1275 


1 


0.010 


0.020 


_ 


61 


126 


202-001+008 


3056 


15 


255.6 


4315.3 


107.0 


-19.25 


3048 


13 


0.130 


0.460 


0.28 


94 


158 


230+027+006 


1830 


10 


215.4 


2263.4 


54.6 


-19.40 


1632 


9 


0.110 


0.399 


0.28 


336 


109 


172+054+007 


1005 


6 


170.0 


1003.6 


53.3 


-19.30 


832 


4 


0.081 


0.250 


0.33 


350 


160 


230+008+003 


955 


5 


105.8 


436.3 


23.0 


-17.50 


947 


2 


0.029 


0.045 


0.65 


38 


95 


167+040+007 


586 


5 


224.4 


660.7 


22.5 


-19.40 


507 


2 


0.028 


0.033 


0.84 


64 


164 


250+027+010 


619 


2 


301.7 


1305.4 


55.7 


-19.70 


618 


3 


0.087 


0.230 


0.38 


87 


_ 


215+048+007 


445 


1 


213.4 


477.8 


21.6 


-19.00 


445 


2 


0.036 


0.033 


1.09 


136 


271 


189+017+007 


504 


5 


212.1 


523.2 


20.7 


-19.25 


446 


2 


0.033 


0.017 


0.97 


152 


160 


230+005+010 


423 


2 


301.6 


907.5 


32.8 


-19.75 


422 


3 


0.068 


0.081 


0.84 


189 


_ 


126+017+009 


433 


1 


267.2 


771.0 


43.7 


-19.70 


409 


4 


0.066 


0.195 


0.35 


198 


82 


152-000+009 


473 


1 


284.7 


863.9 


38.7 


-19.80 


454 


4 


0.056 


0.086 


0.65 


223 


111 


187+008+008 


462 


2 


268.3 


703.7 


34.0 


-19.60 


442 


3 


0.054 


0.139 


0.39 


228 


133 


203+059+007 


643 


4 


210.6 


644.0 


31.1 


-19.10 


612 


2 


0.043 


0.050 


0.86 


317 


_ 


156+010+010 


351 


1 


321.6 


846.7 


36.0 


-20.00 


345 


5 


0.077 


0.226 


0.34 


332 


106 


175+005+009 


333 


1 


291.0 


664.3 


27.3 


-19.90 


309 


3 


0.058 


0.083 


0.69 


349 


138 


207+026+006 


893 


5 


188.0 


768.8 


42.6 


-19.15 


703 


5 


0.056 


0.119 


0.47 


351 


138 


207+028+007 


615 


4 


225.4 


689.1 


33.0 


-19.15 


611 


4 


0.058 


0.088 


0.65 


362 


158 


232+029+006 


306 


2 


195.2 


284.5 


15.1 


-18.80 


305 


2 


-0.004 


-0.004 


1.00 


366 


158 


217+020+010 


353 





300.4 


763.4 


31.1 


-19.90 


339 


4 


0.063 


0.165 


0.38 


376 


167 


255+033+008 


437 


2 


258.7 


658.0 


27.9 


-19.50 


415 


3 


0.058 


0.029 


1.03 


474 


76 


133+029+008 


389 


3 


251.2 


612.6 


43.3 


-19.50 


377 


5 


0.069 


0.221 


0.31 


512 


91 


168+002+007 


371 


2 


227.7 


410.7 


26.7 


-19.40 


321 


3 


0.036 


0.091 


0.39 


525 


109 


177+055+005 


438 


4 


154.2 


312.6 


19.3 


-18.50 


409 


3 


0.017 


0.020 


0.84 


530 




192+062+010 


333 





306.8 


790.3 


40.4 


-20.00 


316 


5 


0.077 


0.223 


0.34 


548 


143 


216+016+005 


314 


2 


158.7 


227.1 


14.6 


-18.50 


294 


1 


0.014 


0.002 




549 




214+001+005 


322 


1 


162.5 


225.8 


16.8 


-18.60 


308 


1 


0.015 


0.027 




550 


154 


227+007+004 


459 


3 


135.1 


287.5 


18.2 


-18.10 


447 


1 


0.004 


0.041 




779 




146+054+004 


353 


1 


139.7 


209.5 


20.4 


-18.00 


353 


2 


0.013 


0.014 


0.94 


796 


93 


167+026+003 


369 


1 


102.3 


161.7 


11.9 


-17.70 


336 


1 


-0.003 


0.001 




827 




189+003+008 


405 





254.1 


572.4 


30.0 


-19.50 


384 


4 


0.045 


0.128 


0.35 



Notes. Columns in the Table are as follows: (1): SCljo' the supercluster ID in the LIO catalogue; (2): IDeoi' the supercluster ID in the EOl 
catalogue; (3): the supercluster ID AAA+BBB+ZZZ, where AAA is R.A., +/-BBB is Dec, and CCC is lOOz; (4): the number of galaxies in the 
supercluster, N^^y, (5): the number of groups with at least 30 member galaxies in the supercluster, N^q', (6): the distance of the density maximum, 
^peak' (^)- '■^^ '■"'■^l weighted luminosity of galaxies in the supercluster, Ljqj; (8): the supercluster diameter (the maximum distance between 
galaxies in the supercluster), Diam; (9): the absolute magnitude limit of the volume limited sample, Rjj^i^gy', (10): the number of galaxies in the 
volume limited sample, A'ggjV' (^ maximum value of the fourth Minkowski functional, ( Vs.max (dumpiness), for the supercluster; (12 ~ 14): 

the shapefinders Ki (planarity) and K2 (filamentarity), and the ratio of the shapefinders Kt /Ki for the full supercluster (the ratio of the shapefinders 
is not determined for superclusters with V3 max = 1 for the full range of threshold densities, as explained in text). 



We calculated the smoothed luminosity density field of 
galaxies and determined extended systems of galaxies (super- 
clusters) using this density field. To determine superclusters, we 
created a set of density contours by choosing a series of density 
thresholds. We define connected volumes above a certain den- 
sity threshold as superclusters. Different threshold densities cor- 
respond to different supercluster catalogues. In order to choose 
proper density levels to determine individual superclusters, we 
analysed the density field superclusters at a series of density 
levels. The mean luminosity density of our sample is (mean - 
1.526-10-2 '°'i^'pf;i . We chose the density level D = 5.0 (in the 
units of mean density) to determine individual superclusters. At 
this density level superclusters in the richest chains of super- 
clusters in the volume under study still form separate systems; at 
lower density levels they merge into huge percolating systems. 



At higher threshold density levels superclusters are smaller and 
their number decreases. Details of the calculation of the lumi- 
nosity density field a nd of the supercluster c atalogue are given 
in Appendix lAl and in lLiivamagi et al.l (l2010h . 

Figure |2] shows the richness of superclusters vs. their dis- 
tance. At distances less than 200 /i"' Mpc there are only a few 
superclusters with less than 300 member galaxies. This selec- 
tion effect is owed to the sample geometry - it is a pyramid with 
the top at the location of the observer At the distance of 100 

Mpc the size of the base is 220x140 /z"' Mpc only, and it 
contains only a few superclusters. At distances between 100 and 
200 /i"' Mpc the number of superclusters is small - there is a 
void region between the nearby superclusters and those at dis- 
tances larger than 200 h^^ Mpc (we will describe the large-scale 
distribution of superclusters in Sect. 4). To avoid these selec- 
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10 



100 200 300 400 
Distance fh"^ MdcI 



500 



Fig. 1. Richness of groups vs. their distance in the TIO group 
catalogue. Only the data for groups with at least six member 
galaxies are plotted. This plot shows that at distances larger than 
» 320 /i ' Mpc almost no rich group is visible. 
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1000 



200 
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Fig. 2. Richness of superclusters vs. their distance (in Mpc//i). 
The line marks the value Ngai - 300. 



tion effects and to be able to resolve the details of supercluster's 
density distribution (their morphology), we chose all superclus- 
ters in our distance interval with at least 300 observed member 
galaxies; Fig.[2]shows that they are present at all distances. Data 
of these superclusters are given in Table [T] In this table super- 
clusters are ordered as they are presented in the text below: first 
the data of the superclusters with at least 950 member galaxies, 
and then of the superclusters with less members. Throughout the 
paper we use the supercluster ID numbers from the LIO cata- 
logue {SCliD in Table[TJ- 

To make the calculations of morphology insensitive to selec- 
tion corrections, we work with volume-limited samples and the 
number density of galaxies instead of flux-limited samples and 
luminosity density. We recalculated the density field for each in- 
dividual supercluster with a kernel estimator with a B3 box spline 
as the smoot hing kernel, with the radius of 8 hr^ M pc (for details 
we refer to Saar et al ] l2007t lEinasto et al.ll2007dl) . and volume- 
limited samples of individual superclusters. The absolute mag- 
nitude limits and the numbers of galaxies in the volume-limited 
versions of superclusters are given in Table [1] 



3. Methods 

3.1. Minkowski functionals and shapefinders 

The supercluster geometry (morphology) is defined by its outer 
(limiting) isodensity surface and its enclosed volume. The mor- 
phology of the isodensity contours is (in the sense of global ge- 
ometry) completely characterised by the four Minkowski func- 
tionals Vo - Vs (we give the formulas in Appendix |B]i. For a 
given surface the four Minkowski functionals (from the first 
to the fourth) are proportional to the enclosed volume V, the 
area of the surface S , the integrated mean curvature C, and the 
integrated Gaussian curvature x- The last of them, the fourth 
Minkowski functional V-i, describes the surface topology; it is a 
sum of the number of isolated clumps and the number of void 
bubbles minus the number of tunnels (voids open from both 
sides) in the region (see, e.g. Saar et al. 2007). High values of the 
fourth Minkowski functional Vt, suggest a complicated (clumpy) 
morphology of a supercluster 

For the argument labelling the isodensity surfaces, we use 
the (excluded) mass fraction mf - the ratio of the mass in re- 
gions with lower density than at the surface, to the total mass 
of the supercluster The value mf - corresponds to the whole 
supercluster, and mf = 1 to its highest density peak. 

With the fourth Minkowski functional V-^ we describe the 
dumpiness of the galaxy distribution inside superclusters - the 
fine structure of superclusters. When the density level is higher 
than the value used to determine a supercluster, the isodensity 
surfaces move from the outer regions of a supercluster into its 
central regions (the value of the mass fraction runs from to 1). 
Therefore some galaxies in the outer regions of a supercluster 
do not contribute to the supercluster any more and this changes 
the inner morphology of the supercluster, which is reflected in 
the V-i - mf relation (the number of isolated clumps changes, 
void bubbles, and tunnels may appear inside superclusters). We 
calculate for superclusters for a range of threshold densities 
(mass fractions), starting with the lowest density that determines 
superclusters (mf = 0), up to the peak density in the supercluster 
core [mf = 1). In Section l4!2l we present figures showing the 
fourth Minkowski functional Vj, for the whole threshold density 
interval for each supercluster and give the maximum value of V3 
for each supercluster. 

The first three Minkowski functionals have been used to 
calculate the dimensionless shapefinders Ki (planarity) an d K2 
(filam entarity) (Sahni et al. 1998; Shandari n et fll. 200ii ISaad 
I2009h . The shapefinders are calculated in two steps. At first, 
specific combinations of Minkowski functionals are used to 
calculate the shapefinders, which describe the thickness, the 
width, and the length of a supercluster With the thickness and 
the width we calculate the planarity Ki of a supercluster, and 
with the width and the length we calculate the filamentarity 
K2 of a supercluster We give their formulae in Appendix |B] 
Einasto et al. (2007d) showed that in the {Ki,K2) shapefinder 
plane the morphology of superclusters is described by a char- 
acteristic curve (morphological signature). When the mass frac- 
tion increases, the changes in the morphological signature ac- 
company the changes of the fourth Minkowski functional. As 
the mass fraction increases, the planarity Ki almost does not 
change for a wide range of mass fractions (up to mf w 0.35), 
while the filamentarity K2 increases - higher density regions of 
a supercluster are more filamentary than the whole supercluster 
At higher mass fractions the planarity of superclusters s tarts to 
decrease. In rich superclusters studied previously (Einas to et al.l 
2OO7J) at a mass fraction of about mf - 0.7, the characteristic 
morphology of the supercluster changes rapidly. Both the pla- 
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Fig. 3. The distribution of groups with at least twelve member galaxies in our superclusters in cartesian coordinates (see text), in 
units of /i ' Mpc. To make projection effects less significant we plot in the right panel a slice with y > 140/1 ' Mpc. The filled circles 
denote groups with at least 50 member galaxies, empty circles denote groups with 30-49 member galaxies and crosses denote 
groups with 12-29 member galaxies. The numbers are supercluster ZD's (Table [1] col. 1). Larger ID numbers show superclusters 
with at least 950 member galaxies, smaller ID numbers - superclusters with less than 950 member galaxies. 



narity and filamentarity decrease - this is a crossover from the 
outskirts of the supercluster to the core of the supercluster In 
high-density, clumpy cores of superclusters, where the isoden- 
sity surfaces have a complex shape, the planarity and filamentar- 
ity m ay even become negative (Sheth et al. 2003; Einasto et al. 
|2008|) . In Fig. l9l(Sect. 14.3b we plot the morphological signature 
with symbols of a size proportional to the value of the fourth 
Minkowski functional V3 at a given mass fraction, to show how 
the dumpiness of the supercluster changes together with the 
change in the morphological signature. 

The morphological signature characterises the morphology 
of superclusters in the whole threshold density interval (the 
fine structure of superclusters). The values of the shapefind- 
ers Ki and K2, and their ratio, K1/K2 (the shape parameter) 
for the whole supercluster quantify the overall shape of su- 
perclusters. The ratio of shapefinders has been used to char- 
acterise the shape of the whole superclus t er, for exaniple, b y 
Kolokotroni s et al ] (l2002l) : lBasilakosl(l2003h : ISheth et alJ(l2003l) : 
Costa-Duart e et al.1 (1201 ll) . 

3.2. Multidimensional normal mixture modelling with Mclust 

We employed multidimensional normal mixture modelling to 
search for possible subsets among superclusters according to 
their physical and morphological parameters. To find an optimal 
model for the collection of subsets, the Mclust package for clas- 
sification and clustering was applied. This package searches for 
an optimal model for the clustering of the data among models 
with varying shape, orientation and volume, finds the optimal 
number of concentrations, and the corresponding classification 
(the membership of each concentration). 

The Mclust package gives two statistical measures to esti- 
mate how well the superclusters are divided into the subsets. 
First, Mclust calculates the classification uncertainty of for each 
object in a dataset. This parameter is defined by the probabilities 
for each object to belong to a particular subset and is calculated 
as 1 . minus the highest probability of a supercluster to belong 
to a given subset. The classification uncertainty can be used as a 
statistical estimate of how well objects are assigned to the sub- 
sets. 



To measure how well the subsets are determined and to 
find the best model for a given dataset, Mclust uses the Bayes 
Information Criterion (BIC), which is based on the maximized 
log-likelihood for the model, the number of variables and the 
number of mixture components. The model with the lowest value 
of the BIC among all models c alculated by Mclust is con sidered 
the best. For details we refer to fealev & Raftervl(l2006l) . Below 
we calculate both the uncertainities of the classification of su- 
perclusters in the best model determined using Mclust, and the 
values of the BIC for different classifications of superclusters as 
found by Mclust. 



4. Results 

4. 1 . Large-scale distribution of superclusters 

We start with the identification of Abell clusters among groups 
with at least 30 member galaxies in our superclusters. With these 
data we identify our superclusters with those determined earlier 
on the basis of Abell clusters (EOl). This will help to analyse the 
large-scale distribution of superclusters and to compare it with 
earlier studies. 

We present a list of the Abell clusters in superclusters in 
Table |2] Here the X-ray clusters are also marked; about 1 /3 of 
the Abell clusters in Table|2]are X-ray sources. These data were 
used to compare superclusters with those determined in EOl . We 
give the EOl's ID number if there is at least one Abell cluster in 
common between EOl and the present supercluster sample. A 
word of caution is needed - a common cluster does not always 
mean that superclusters can be fully identified with each other 
A number of superclusters from EOl are split between several 
superclusters in our present catalogue. In these cases the iden- 
tification of superclusters with the superclusters found on the 
basis of the Abell clusters is complicated and has to be taken as 
a suggestion only. 

Table [T] and table |2] show that in our sample there are three 
superclusters without any group/cluster with at least 30 member 
galaxies, and seven superclusters that do not contain Abell clus- 
ters. Among the superclusters eight have one group/cluster with 
at least 30 member galaxies. All these systems are poor, compa- 
rable with the Local Supercluster, which has only one rich clus- 
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ter. Two of the superclusters contain at least 10 groups/clusters 
with at least 30 member galaxies; these are the richest and most 
luminous superclusters in the sample: SCI 061 and the Corona 
Borealis. 

Table 2. Abell clusters in superclusters. 



(1) 


(2) 


(3) 


S Clio 


IDeoi 


Abell ID 


1 


162 


2142x, 2149x 


10 


160 


2152 


11 


154 


2040, 2028 


24 


111 


1424, 1516 


38 


95 


1173x, 1187 


55 


111 


1358 


60 


160 


2197 


61 


126, 136 


1620, 1650x, 1658x, 1663x, 1692, 






1750x, 1773x, 1780, 1809x 


64 


164 


2223, 2244 


87 




1904 


94 


158 


2067, 2065x, 2089 


136 


271 


1569 


152 


160 


2048, 2055x 


198 


82 


933 


223 


111 


1541, 1552 


228 


133 


1767x 


332 


106 


1346 


336 


109 


1279, 1436 


349 


138 


1795x, 1827, 1831x 


350 


160 


2052x, 2063x 


351 


138 


1775x, 1800x, 1831x 


362 


158 


2073, 2079, 2092 


376 


167 


2249x 


474 


76 


699 


512 


91 


1205x, 1238 


525 


109 


1291x, 1377 


548 


143 


1913 


550 


154 


2028, 2055 


796 


93 


1185x 



Notes. Columns in the table are as follows: (1): 5 Cho'- the supercluster 
ID in the LIO catalogue; (2): /D^oi : the supercluster ID in the EOl cata- 
logue ; (3): the Abell ID. x denotes X-ray clusters ("EOl. lBohringer et al.l 
l2004h . 



We show the large-scale distribution of rich clusters in su- 
perclusters in Fig.|3]in cartesian coordinat es. These coordinates 
are defined as iu iPark et al.. (.2007.) : .Liivamagi et al.l (i2010i) : 

X - -d sin A, 

y - dcosAcosT], (1) 
z = li cos /I sin 7/, 

where d is the comoving distance, and A and rj are the SDSS 
survey coordinates. To complement this figure, we present in 
Appendix ICl a 3D version of Fig. [3] and the 3D distributions of 
groups in superclusters with the right ascensions, declinations, 
and distances of groups. 

In Fig.[3]the superclusters SCI 060 and SCI 350 (Table[D are 
seen close to us. They belong to the Hercules supercluster, which 
is split into several superclusters in our sample. These are the 
nearby rich systems seen in Fig.|2] A chain of poor superclusters 
connects the Hercules supercluster with rich superclusters at a 
distance of about 200 Mpc the superclusters SCI 349 and 
SCI 351 (the Bootes supercluster) among them. We mentioned 
in Sect. 2 that at small distances the size of the sample cross- 
section is only 220x140 /i"' Mpc and these superclusters may be 



broken up by the sample borders. Therefore the data on nearby 
superclusters are less reliable than the data on the more distant 
ones. 

Rich superclusters at distances of about 210-260 /i"' Mpc 
form three chains, separated by voids. A 3D figure on our web 
pages shows that actually only one of these supercluster systems 
ia a clear chain-like system. This is the Sloan Grea t Wall (SGW), 



2004: Gottetal. 20051 Nicholetal, 


200d: Einasto et al. 


|20ld; 


Luparello et al.ll201 ll iPimbblet et al. 


201lHEinasto et al.1 


201 ih. 



The SGW consists of several superclusters of galaxies. The rich- 
est of them are the superclusters SCI 061 and SCI 024. The other 
two supercluster chains are much poorer and cannot really be 
called chains. One of them is separated from the SGW by a void; 
the Bootes supercluster is the richest supercluster in this system. 
The richest supercluster in the third system of superclusters is 
the Ursa Majoris supercluster (SCI 336). 

The very rich Corona Borealis supercluster is located at 
the joint of these systems. This supercluster is a member of a 
huge system of rich superclusters located at the ri ght angle with 
respect to the Local Supercluster, described by lEinasto et al.l 
(Il997bh as the dominant supercluster plane. 

At high positive values of the x coordinate there are no rich 
clusters, and the superclusters in this region are also poor. There 
are some poor superclusters farther away, perhaps connecting 
superclusters in our sample volume with more distant superclus- 
ters. Thus the large-scale distribution of the superclusters is very 
inhomogeneous, as noted also in LIO. 

To make projection effects less significant, we do not show 
the chains of nearby superclusters in the right panel of Fig. [3] 
The superclusters SCI 548, SCI 549, SCI 550, SCI 350, and 
SCI 060 are superimposed on the superclusters SCI 351 and 
SCI 354, and the superclusters SCI 779 and SCI 796 on the su- 
perclusters SCI 038, SCI 336, and SCI 525. 

4.2. Morphology of superclusters 

The results on the morphology of superclusters are summarised 
in Table [1] where we list the following morphological character- 
istics for each supercluster: the maximum value of the fourth 
Minkowski functional Vs.max (dumpiness), the values of the 
shapefinders Ki (planarity) and K2 (filamentarity) for the whole 
supercluster, and the ratio of the shapefinders K1/K2 for the 
whole supercluster This ratio is not given for the superclusters 
for which V3 = I over the whole mass fraction interval, because 
the K2 may become very small, making the ratio K1/K2 noisy. 

We began the analysis of the structure of superclusters by 
searching for possible subsets defined by their physical and 
morphological characteristics. For this purpose we applied the 
Mc/M.sf package for classification and clustering, decribed shortly 
in Sect. 13.21 Initial data for Mclust was the number of groups 
with at least 30 member galaxies, the total weighted luminosity 
of galaxies in a supercluster, its diameter, and the morphological 
parameters given in Table [1] The results of this analysis show 
that superclusters can be divided into two main sets. The first set 
consists of superclusters with shape parameter K1/K2 < 0.6, and 
second set of those with shape parameter K1/K2 > 0.6 - these 
superclusters are less elongated than the superclusters in the first 
set. According to Mclust, the richest supercluster in our sam- 
ple, the supercluster SCI 061, forms a separate subgroup owing 
to its very high luminosity. For simplicity, we include this su- 
percluster in the following analysis in the set of more elongated 
superclusers. 
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We estimated how well superclusters are assigned to the two 
different sets using the uncertainity of classification calculated 
by Mclust (see Sect. l3.2l i. For our superclusters, the mean uncer- 
tainty of the classification is 1 .9 • 1 0"^, showing that superclusters 
are well classified. To measure how well the sets are determined 
and to find the best model of a given dataset we use the BIC 
values given by Mclust. For our sets of superclusters the low- 
est value of the BIC corresponds to the division of superclusters 
into two main sets, with SCI 061 forming a separate class. The 
value of the BIC for the one-component model for superclusters 
is higher, showing that this model is less likely for our superclus- 
ters. 

Table 3. Parameters of superclusters in sets according to their 
morphology. 



(1) 


(2) 


(3) 


(4) 




K1/K2 


< 0.6 


K, IK2 > 0.6 




with SCI 061 


without SCI 061 






16 


15 


20 


Distance 


254 ± 64 


254 ± 66 


200 ± 43 


^tot 


820 + 316 


790 ± 274 


525 ± 122 




530 + 207 


442 ± 176 


445 ± 129 


^30 


2.5 ± 0.92 


2.0 ± 0.74 


2.0 ± 0.63 


Diameter 


44 ± 12 


43 ± 11 


22 ± 5 


Vi 


4.5 ± 1.25 


4.0 ± 1.15 


2.0 ± 0.48 


Ki 


0.08+ 0.02 


0.08 ± 0.02 


0.03 ± 0.007 


K2 


0.20± 0.05 


0.200 ± 0.05 


0.031 ± 0.008 



Notes. Columns in the Table are as follows: 1: The parameter: A^s^, 
the number of superclusters; the median distance of superclusters, in 
Mpc//z; the median total weighted luminosity of galaxies in the super- 
clusters, i]jjg£j, in IO'^/i^-Lq; the median number of galaxies in the 
superclusters, N^^ \ the median number of groups with at least 30 mem- 
ber galaxies in the superclusters, A'3o; the median diameter of the su- 
perclusters, in Mpc//7; the median value of the fourth Minkowski func- 
tional, V3; the median values of the planarity Ki and filamentarity K2 
of the superclusters. 2-3: The median values and Icr errors of the cor- 
responding parameter for the first set of superclusters with the ratio of 
the shapefinders KJK2 < 0.6; 4: The median values and Icr errors of 
the corresponding parameter for the second set of superclusters with the 
ratio of the shapefinders KJK2 > 0.6. 



Fig. IHpresents the shapefinders Ki (planarity) and K2 (fila- 
mentarity) in the shapefinder's plane for our supercluster sam- 
ple. In this figure circles correspond to those superclusters for 
which the ratio of the shapefinders K\ /K2 < 0.6, and squares to 
the superclusters with Ki /K2 > 0.6, i.e. to more elongated and 
to less elongated superclusters, respectively. The symbol sizes 
are proportional to the maximum value of the fourth Minkowski 
functional VB^max (dumpiness). This figure reflects both the outer 
shape and the inner dumpiness (fine structure) of superclusters 
and summarises the morphological information about superclus- 
ters. 

Table [1] and Fig. |4] demonstrate that almost all superclus- 
ters in our sample are elongated; they have large filamentari- 
ties with larger range of values than planarities. Two superclus- 
ters with the largest filamentarities and with the ratio of the 
shapefinders as small as K1/K2 - 0.28 are the most extreme 
cases of filamentary systems in our sample. These two super- 
clusters are the richest in our sample: SCI 061 and the Corona 
Borealis. We do not have extremely planar superclusters in our 
sample (K1/K2 » 1.0). Three superclusters in our sample have 
Ki IK2 ~ 1 .0 (they are spherical), and all these have a small max- 
imum dumpiness (Vs^max ^ 3). More elongated superclusters 



CM 




0.15 



K1 



Fig. 4. Shapefinders K\ (planarity) and K2 (filamentarity) for the 
superclusters with their dumpiness. The symbol sizes are pro- 
portional to the fourth Minkowski functional Va.max- Circles de- 
note superclusters with the shapefinders ratio K1/K2 < 0.6 and 
squares denote the superclusters with the ratio K1/K2 > 0.6. 
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Fig. 5. Cumulative distributions of the physical characteristics 
of superclusters for the two sets of superclusters divided by the 
value of the shape parameter. From left to right: the distributions 
of luminosities, numbers of galaxies, numbers of rich groups 
with at least 30 member galaxies, and diameters of superclusters. 
The black lines correspond to the first set of superclusters with 
K1/K2 < 0.6, including the supercluster SCI 061, the grey lines 
to the first set of superclusters without SCI 061, and the dashed 
lines to the second set of superclusters with K1/K2 > 0.6. 



typically have higher values of dumpiness V3 than less elon- 
gated ones - they have a more complicated morphology (see 
Table[3]l. 

We present the median values of supercluster parameters in 
these sets in Table [3] In Fig.|5]we plot the cumulative distribu- 
tions of the values of the physical characteristics of superclusters 
from the two sets, and in Fig. |6]the cumulative distributions of 
the morphological parameters. The scatter of the parameters is 
large, but the superclusters in the first set with the shape parame- 
ter K1/K2 < 0.6 (with and without the supercluster SCI 061) are 
richer, more luminous, and have larger diameters than those in 
the second set with shape parameter K1/K2 > 0.6. Superclusters 
from the first set also have higher maximum values of the fourth 
Minkowski functional V3, and higher values of planarities Ki 
and filamentarities K2. In Fig. [T] we show another presentation of 
these results - the shapefinder's plane for superclusters of differ- 
ent total luminosity and the richness (the number of galaxies in 
a supercluster) as coded in the symbol sizes explained in figure 



7 



M. Einasto et al.: Morphology of superclusters 





Fig. 6. Cumulative distributions of morphological parameters of 
superclusters for two sets of superclusters divided by the value 
of the shape parameter From left to right: the distributions of 
the fourth Minkowski functional Vj, of planarities Ki, and fil- 
amentarities K2 of superclusters. The black lines correspond to 
the first set of superclusters with K1/K2 < 0.6, including the 
supercluster SCI 061, the grey lines to the first set of superclus- 
ters without SCI 061, and the dashed lines to the second set of 
superclusters with K1/K2 > 0.6. 
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Fig. 7. Shapefinders Ki (planarity) and K2 (filamentarity), with 
the luminosities and richnesses of the superclusters. In the left 
panel symbol sizes are proportional to the total weighted lumi- 
nosity of galaxies in the superclusters, and in the right panel to 
the number of galaxies in the superclusters. Circles and squares 
correspond to two sets of superclusters as in Fig. 2] 



captions. The two most elongated superclusters are the richest 
and the most luminous. 

Figure |8] and Figure |5] show that more elongated superclus- 
ters also have larger diameters than less elongated ones, imply- 
ing that the systems with larger diameters are not planar struc- 
tures. Also, there are no compact, planar, and very luminous su- 
perclusters. 

Table [T] shows that all superclusters with the ratio of 
shapefinders with K1/K2 > 0.6 are relatively nearby, only four 
of them lie at distances larger than 250 /i"' Mpc. The closest 
of them is located at about 90 /i ' Mpc, and their median dis- 
tance is about 200 /i"' Mpc. As mentioned above, owing to the 
sample geometry, the nearby superclusters may not be fully in- 
cluded in the sample volume and their small dumpiness may be 
due to this selection effect. Some nearby poor superclusters are 
located in low-density filaments between us and more distant su- 
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Fig. 8. Diameters of superclusters (in Mpc//i) vs. their shape pa- 
rameter Symbol sizes are proportional to the value of the fourth 
Minkowski functional of a supercluster Circles and squares 
correspond to the two sets of superclusters as in Fig. 



perclusters, a nd their shape and sm all dumpiness may be real (as 
noted also in lEinasto et aLlllQgfb ). The closest more elongated 
superclusters with the shape parameter K1/K2 < 0.6 are about 
170 /i"' Mpc away from us, their mean and median distances al- 
most coincide and are about 254 /z ' Mpc (this is approximately 
the distance to the rich superclusters in the SGW). 

4.3. Notes on individual superclusters 

Below, we give a short description of individual superclusters in 
our sample. In the figures of this section and in Appendix ID] we 
show for each supercluster (except for those for which Vy, = I 
over the whole mass fraction interval) the sky distribution of su- 
percluster members, the values of the fourth Minkowski func- 
tional V3 vs. the mass fraction mf and the morphological sig- 
nature for each supercluster Panels in these figures are as fol- 
lows. The left panels show the sky distributions of galaxies in 
superclusters and the location of rich groups with at least 30 
member galaxies. The middle panels show the dumpiness V3 
vs. the mass fraction mf for a whole mass fraction interval from 
to 1, and the right panels the shapefinder's {K[,K2) curve (the 
morphological signature) for a supercluster The mass fraction 
increases anti-clockwise along the curves. In these panels the 
value of mass fraction mf — 0.7 is marked - at this value the 
morphological signature of rich superclusters changes. We will 
classify our superclusters as spiders, multispiders, filaments, and 
multibranching filaments on the basis of their morphological in- 
formation and visual appearance. Often superclusters are of in- 
termediate type between these main types, hence for some su- 
perclusters our classification is a suggestion only. 

We begin with the two richest superclusters from our sample 
(Fig.|9]l. For these two superclusters the morphological signature 
(Fig S right panels) is plotted with symbols of the size propor- 
tional to the value of the fourth Minkowski functional Vt, at a 
given mass fraction, to show how the morphological signature 
changes together with the changes in dumpiness, as described 
in Sect. 13.11 Here we do not mark the value of mass fraction 
mf - 0.7, for clarity. 

The supercluster SCI 061 at a distance of 25 6 Mpc is the 
richest member of the SGW (Einasto et al. 201 l: iLuparello et al.l 
201 li) . This supercluster contains nine Abell clusters, the largest 
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Fig. 9. Left panels: the distribution of galaxies (grey dots) in the sky, for the two superclusters, SCI 061 (upper row) and SCI 094 
(lower row). Circles mark the location of groups with at least 30 member galaxies, and the size of a circle is proportional to the size 
of a group in the sky. The numbers show Abell clusters. The middle panels show the fourth Minkowski functional V^, and the right 
panels the shapefinders Ki (planarity) and K2 (filamentarity) for a supercluster. The morphological signature in the Ki - K2 plane is 
parametrically defined as Kiim f) and K2{m f). The sizes of open circles are proportional to the value of at a given mass fraction 
mf. They show the change of the dumpiness with the mass fraction together with the changes in the morphological signature. The 
Abell clusters in the supercluster SCI 061 (upper row) are 1 - A1620, 2 - A1650, 3 - A1658, 4 - 1663, 5 - 1682, 6 - 1750, 7 - 
1773, 8 - 1780, and 9 - 1809. The Abell clusters in the supercluster SCI 094 (lower row) are 1 - A2067, 2 - A2065, and 3 - A2089. 



number in our sample. Five of these are also X-ray clus- 
ters (Bohringer et al. 2004). The richest of them, A1750, is a 
merging X-ray cluster (Belsole et al. 2004). The morphology of 
SCI 061 resembles a multibranching filament with the maximum 
value of the fourth Minkowski functional V3.max = 13, and the 
ratio of the shapefinders K1/K2 = 0.28, one of the lowest in our 
catalogue (Fig.|9] upper row, and Table [1}. 

The supercluster SCI 094 (the Corona Borealis supercluster) 
at a distance of 215 h'^ Mpc is the second in richness among 
our sample. This system contains three Abell clusters (A2067, 
A2065, and A2089), and is a mem ber of the dominant super- 
cluster plane dEinasto et al.lll997bl) . The distribution of galax- 
ies in the sky in SCI 094 is plotted in Fig. |9] (lower row, left 
panel). The maximum value of the fourth Minkowski functional 
of the Corona Borealis supercluster Vs max = 10, and the ra- 
tio of the shapefinders K1/K2 - 0.28 (Fig. |9] the middle and 
right panels of the lower row, and Table [1}. Morphologically 
SCI 094 is a multispider with a number of clusters connected 
by low density filaments, with an overall very elongated shape 
that resembles a horse-shoe, with the merging X-ray cluster 
A2065 at the top ( Chatzikos et al. 2006). SCI 094 has been 
studied by ISmall et al] (1 1998b . who found that the core of this 
system probably has started to collapse. Numerical simula- 



tions show that such collapsing cores in superclusters are rare 
(Gramann & Suhhonenko 2002). Luparello et al. (2011) pro- 
posed that SCI 094 may merge with several surrounding su- 
perclusters in the future. In last years interest in the Corona 
Borealis supercluster region has grown because of the discov- 
ery of the CMB cold spot in it's direction. This may be partly 
caused by the warm-hot diff'use gas in the supercluster fila- 
ments b etween the clusters, or by some undiscovered distant 
clus ter (IGenov a-Santos eTal? '2008'; 'Padil la-Torres et al.l 120091: 
Gen ova-Santo s et al. 2010; Padilla-Torres et al.l2010l and refer- 
ences therein). The poor superclusters SCI 362 and SCI 366 are 
also members of the Corona Borealis supercluster. The morphol- 
ogy of SCI 362 resembles a simple spider with the ratio of the 
shapefinders K1/K2 - 1.0, while that of SCI 366 resembles a 
multibranching filament with the maximum value of the fourth 
Minkowski functional Vj, max - 4, and the ratio of the shapefind- 
ers K1/K2 =0.38 (Fig. EIJ- 

The supercluster SCI 001 contains two rich Abell clusters, 
A2142 and A2149, both of them are X-ray sources (Table |2). 
Chandra observations have revealed that A2142 is probably still 
merging (Markevitch et al. 2000). SCI 001 is located in a region 
with a dense concentration of superclusters, close to SCI Oil and 
SCI 094. All these systems form a part of the dominant super- 



9 



M. Einasto et al.: Morphology of superclusters 



cluster plane. At the location of SCI 001 the luminosity density 
is the highest in the whole SDSS survey; this is probably at least 
partly because of the rich X-ray cluster A2142. The morphol- 
ogy of SCI 001 resembles a filament where clusters are located 
almost along a straight line (the 3D model on our web pages 
shows this best). The maximum value of the fourth Minkowski 
functional for SCI 001 Vj max - 2 and the ratio of the shapefind- 
ers K1/K2 = 0.48 (Fig.[TUiand TableQ. 

The supercluster SCI Oil at a distance of 234 /i"' Mpc con- 
tains two Abell clusters, A2040 and A2028, which are mem- 
bers of different superclusters in EOl . The maximum value of the 
fourth Minkowski functional Vs ^ax = 3 (Fig. [TOl middle row), 
the ratio of the shapefinders K1/K2 = 0.92. The morphology of 
SCI Oil resembles a sparse multispider or multibranching fila- 
ment with a quite uniform density (as suggested by the low val- 
ues of the fourth Minkowski functional V3 for a wide mass frac- 
tion interval). SCI Oil belongs to the same supercluster complex 
as SCI 001 (Fig.EJ. 

The supercluster SCI 024 at a distance of 230 Mpc is 
the second richest member of the SGW. SCI 024 contains two 
Abell clusters, A1424 and A1516. Its morphology resembles a 
multispider with the maximum value of the fourth Minkowski 
functional y3,,nax = 5 and the ratio of the shapefinders K1/K2 - 
0.48 (Fig. [To] lower row). The values of the fourth Minkowski 
functional of SCI 024 at small mass fractions (mf w 0.25) 

= suggest th at SCI 024 has low-density tunnels inside 
jEinasto et al.ll20lTl) . The member of the supercluster SCI 111 
in EOl to which SCI 024 also belong, the poor supercluster 
SCI 223, can be described as a multispider with two concen- 
trations -the fourth Minkowski functional has a value of 2 for 
a wide mass fraction interval (Fig. ID.U . SCI 223 is separated 
from SCI 024 and SCI 055 by a small void, showing that clus- 
ters gathered together into one supercluster in the EOl catalogue 
sometimes do not belong to the same supercluster, when systems 
are determined using data on galaxies. 

The supercluster SCI 055 at a distance of 242 h^^ Mpc is 
separated from the SGW by a void and is connected to it by a 
filament of galaxies (this is the reason why SCI 024 and SCI 055 
both belong to the supercluster SCI 111 in the EOl catalogue). 
SCI 055 contains one Abell cluster, A1358. The morphology of 
SCI 055 resembles a multibranching filament or an elongated 
multispider with the maximum value of the fourth Minkowski 
functional V3 max = 5 and the ratio of the shapefinders K1/K2 - 
0.43 (Fig.inii. 

The supercluster SCI 336 contains two Abell clusters, A 1279 
and A 1436, which are members of the Ursa Majoris superclus- 
ter (Fig. [1X1 middle row). SCI 336 belongs to a chain of super- 
clusters that is separated by a void from SCI 349 and SCI 351 



(the Bootes supercluste r) (see also Kopvlova & Kopvlovl [ 20091 
and references therein). Luparello et al.l ( 201 ll) found that sev- 



eral filamentary systems may be associated with this system. 
According to our calculations, the morphology of the Ursa 
Majoris supercluster resembles a sparse multibranching filament 
with the maximum value of the fourth Minkowski functional 
V3,max - 4, and the ratio of the shapefinders K1/K2 - 0.33 
(Fig. fTTT l. The poor supercluster SCI 525 is also a member of 
the Ursa Majoris supercluster, morphologically SCI 525 can be 
described as a spider. 

The superclusters SCI 010, SCI 060, and SCI 350 belong to 
the rich Hercules supercluster that is split between several sys- 
tems in our present catalogue. They are located at a distance of 
about 100 Mpc and contain clusters that are exceptionally 
rich in the TIO catalogue, owing to their small distance from us 
(A2152 in SCI 010 and A2197 in SCI 060). This is the reason 



why the value of the fourth Minkowski functional for SCI 010 
and SCI 060 V3 = 1 over the whole mass fraction interval, show- 
ing that they contain only one high density clump. SCI 350 con- 
tains five rich clusters, three of them correspond to the Abell 
cluster A2052. This suggests that A2052 has a substructure, with 
different components corresponding to different groups and clus- 
ters in the TIO catalogue. The maximum value of the fourth 
Minkowski functional for SCI 350, Vs.max = 2 (Fig. (TT] lower 
row), and the ratio of the shapefinders K1/K2 - 0.65. At the 
farther end SCI 350 chain joins the dominant supercluster plane 
(Fig.EJ. 

We summarise the results of the morphological analysis of 
the superclusters in Table |4] For clarity we present the figures 
of the distribution of galaxies and rich groups/clusters in the sky 
for superclusters with less than 950 members galaxies, as well as 
their fourth Minkowski functionals and the morphological signa- 
tures in Appendix Id] 

In Table |4] we list for every supercluster their distance, the 
numbers of galaxies, the number of rich groups with at least 30 
member galaxies, and the number of Abell clusters among them. 
From morphological parameters we list the maximum value of 
the fourth Minkowski functional and the number that shows 
whether the supercluster belongs to the set 1 (more elongated 
superclusters) or to the set 2 (less elongated superclusters). We 
also give morphological descriptions of superclusters and notes. 

Table |4] shows that superclusters with a smaller number of 
galaxies also contain, as expected, a smaller number of rich 
groups and Abell clusters. Among them less elongated super- 
clusters dominate over more elongated superclusters - there are 
16 systems from set 2 and 10 from set 1 among them. The 
morphology of 14 of them can be described as simple spider 
or simple filament, and 12 of them are either multispiders or 
multibranching filaments. In contrast, of 10 superclusters with 
at least 950 member galaxies 7 can be described as multispiders 
of multibranching filaments, two of them are simple spiders and 
one is a simple filament. 

5. Discussion 

5.1. Selection effects 

The main selection effect in our study comes from the use of the 
flux-limited sample of galaxies to determine the luminosity den- 
sity field and superclusters. To keep the luminosity-dependent 
selection effects as small as possible, we used data on galaxies 
and galaxy systems for a distance interval 90-320 /i"' Mpc. In 
this interval these effects are the smallest (we refer to TIO for 
details). We calculate Minkowski functionals of individual su- 
perclusters from volume limited samples. This approach makes 
the calculations of morphology insensitive to luminosity depen- 
dent selection effects. 

Another selection effect comes from the choice of the den- 
sity level used to determine superclusters. At the density level 
applied in the present paper [D = 5.0), rich superclusters do 
not percolate yet. For example, in the SGW we see several indi- 
vidual rich superclusters. At a lower threshold density these su- 
perclusters join into a huge system. In the same time the Corona 
Borealis supercluster is split into several superclusters in our cat- 
alogue at the density level D = 5.0 . Thus there is no unique den- 
sity level, which would be the best choice for all superclusters. 
However, LIO show that superclusters are fairly well-defined 
systems and do not change much when changing the density 
level, if only this change is small and does not break them up into 
small systems or does not join them into huge percolating sys- 
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R.A. mf K1 



Fig. 10. Panels as in Fig.|9] Filled circles in the right panel mark the value of the mass fraction mf - 0.7. Upper row: the supercluster 
SCI 001. The Abell clusters are 1 - A2142 and 2 - A2149. Middle row: the supercluster SCI Oil. The Abell clusters are 1 - A2040 
and 2 - A2028. Lower row: the supercluster SCI 024. The Abell clusters are 1 - A1424 and 2 - A1516. 



tems. When we move towards higher density levels, some galax- 
ies and galaxy systems in the lower density outskirts of a super- 
cluster do not belong to the supercluster any more; in that case 
the dumpiness of the supercluster at low mass fractions may de- 
crease. Moving towards lower density levels, new galaxies and 
galaxy systems in the outskirts of a supercluster become super- 
cluster members. This may increase the value of the dumpiness 
and change the morphological signature of a supercluster at low 
mass fractions. The possible change in morphology depends on 
the number of galaxies and on the richness of groups removed 
from or added to the supercluster and is individual for each su- 
percluster. 

At small distances our sample volume is small, which in- 
troduces another selection effect - nearby rich superclusters are 
split between several superclusters, some parts of these super- 
clusters may be located outside of the sample volume. This 
makes the data about nearby superclusters less reliable. 



The identification of Abell clusters may also be affected 
by selection effects. It is well known that some Abell clusters 
consist of several line-of-sight components. A good example 
of such a cluster is th e cluster Abell 1386, recently studied by 
Pim bblet et al ] (1201 Ih . Our search discarded Abell clusters af- 
fected by this projection effect. 

5.2. Comparison with otiner studies 

Recently ICosta-Duarte et al" I (l20Tll) used volume-limited sam- 
ples of galaxies from the SDSS DR7 to extract superclusters 
of galaxies and to study the morphology of whole superclus- 
ters with the shape parameter (the ratio of the shapefinders 
K[/K2). To determine superclusters the authors calculated the 
density field with an Epanechnikov kernel and f ound systems 
of galax ies with at least 10 member galaxies. In lEinasto et al.l 
(l2007dh . we compared the Epanechnikov and B3 box spline ker- 
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Fig. 11. Panels are the same as in Fig. |9] Upper row: the supercluster SCI 055. The Abell cluster is A1358. Middle row: the 
supercluster SCI 336. The Abell clusters are 1 - A1279, 2 - A1436. Lower row: the supercluster SCI 350. The Abell clusters are 1 
-A2052,2-A2063. 



nels and found that both kernels are good to describe the over- 
all shape of superclusters, while the Bj, box spline kernel bet- 
ter resolves the inner structure of superclusters. This is the rea- 
son why we used this kernel in the present study. The most im- 
portant difference between our studies is that we used the data 
abo ut the richer superclusters with at least 300 member galax- 
ies. [Costo^Du^teeFal] (|201ll) showed that there are both planar 
and filament-like superclusters (pancakes and filaments) among 
the superclusters of thek sample, while in our sample there are 
almost no superclus ters for which planarity is larger than fila- 
mentarity. However, ICosta-Duarte et all (1201 ll) found that very 
rich and luminous superclusters tend to be filaments, as we also 
found in our study. 

The overall shapes of superclusters, described by the shape 
parameters or approximated by triaxial ellipses, have been 
analysed in [Jaaniste et al. (1998); Kolokotronis et al. (200^ 
iBasilakosI (l2003i) : £inasto et aL (i2007a) : .Costa-Duarte et aP 



(1201 lh : lLuparello et al.1 (l201 ih . These studies showed that elon- 
gated, prolate structures dominate among superclusters, as 
we also found in our study. In the analysis of the geome- 
try of the structures around the galaxy clusters from simu- 
lations iNoh & CohnI (1201 II) showed that these structures tend 
to lie on planes. Similarly, in our study and in the study 
by Costa-Duarte e t al.l (1201 ih poorer systems are more planar. 
Sheth et alj (i2003i) used the shapefinders plane Ki -K2 (the mor- 
phological signature in our study) to study the morphology of 
superclusters in simulations. They found that richer and more 
massive superclusters tend to be more filamentary and have a 
more complicated inner structure (high values of the genus in 
their study). 

Future evolution of the structure in the Universe has 
been addressed in simulations by several authors, we re- 
fer to lAr aya-Melo et al. (2009) for a review and references. 
lArava-Melo et al.i (i2009>) studied the evolution of the shape and 
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Table 4. Data of the morphology of superclusters. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


5 CliD 


"^peak 
Mpc//j 


^gal 


^30 


^Abell 


Vs.max 


Set 


Morph. 


Notes 


1 


264.5 


1038 


4 


2 


2 


1 


F 




10 


111.9 


1463 


2 


1 


1 


2 


S 




11 


235.0 


1222 


5 


2 


3 


2 


mS(mF) 




24 


230.4 


1469 


4 


2 


6 


1 


mS 


SGW 


55 


242.0 


1306 


4 


1 


5 


1 


mS(mF) 




60 


92.0 


1335 


3 


1 


1 


2 


S 




61 


255.6 


3056 


15 


9 


13 


1 


mF 


SGW 


94 


215.4 


1830 


10 


3 


9 


1 


mS 




336 


170.0 


1005 


6 


2 


4 


1 


mF 




350 


105.8 


955 


5 


2 


2 


2 


mS 




38 


224.4 


586 


5 


2 


2 


2 


S 




64 


301.7 


619 


2 


2 


3 


1 


F 




87 


213.4 


445 


1 


1 


2 


2 


S 




136 


212.1 


504 


5 


1 


2 


2 


S 




152 


301.6 


423 


2 


2 


3 


2 


F 




189 


267.2 


433 


1 





4 


1 


F 




198 


284.7 


473 


1 


1 


4 


2 


mF 




223 


268.3 


462 


2 


2 


3 


1 


mS 




228 


210.6 


643 


4 


1 


2 


2 


F 




317 


321.6 


351 


1 





5 


1 


mS 


T 


332 


291.0 


333 


1 


1 


3 


2 


mS 




349 


188.0 


893 


5 


3 


5 


1 


mS 


F(Lll) 


351 


225.4 


615 


4 


3 


4 


2 


mF 




362 


195.2 


306 


2 


3 


2 


2 


mS 




366 


300.4 


353 








4 


1 


mF 




376 


258.7 


437 


2 


1 


3 


2 


F 


T 


474 


251.2 


389 


3 


1 


5 


1 


F 




512 


227.7 


371 


2 


2 


3 


1 


mF 




525 


154.2 


438 


4 


2 


3 


2 


mS 




530 


306.8 


333 








5 


1 


mS 


T 


548 


158.7 


314 


2 


1 


1 


2 


S 


Her-DSP 


549 


162.5 


322 


1 





1 


2 


S 


Her-DSP 


550 


135.1 


459 


3 


2 


1 


2 


S 


Her-DSP 


779 


139.7 


353 


1 





2 


2 


mS 


SGW chain 


796 


102.3 


369 


1 


1 


1 


2 


S 


SGW chain 


827 


254.1 


405 








4 


1 


F 





Notes. Columns in the Table are as follows: (1): SCljo' the supercluster ID in the LIO catalogue; (2): the distance of the density maximum, d^^g^; 
(3): the number of galaxies in the supercluster, A'ggj; (4): the number of groups with at least 30 member galaxies in the supercluster, N^q; (5): 
the number of Abell clusters among groups with at least 30 member galaxies in the supercluster, A'^bell' maximum value of the fourth 

Minkowski functional, (Vsi^^x (dumpiness), for the supercluster; (7): classification: set - 1 for superclusters with the ratio of the shapefinders 
Ki/K2 < 0.6 and 2 for superclusters with the ratio of the shapefinders KijKi > 0.6); (8): morphology: S denotes spiders, mS - multispiders, F - 
filaments and mF - m ultibranching filaments; (9): notes. SGW denotes members of the Sloan Great Wall; F(Lll) denotes a filamentary system 
according to jLupar ello et al. 2011); T denotes supercl usters for which at a certain mass fraction V3 = 0, suggesting that these superclusters 
have tunnels through them (App. IB] .Einasto et al.ll201 lb . Her-DSP denotes the superclusters located in a chain of superclusters that connects the 
Hercules supercluster with the dominant supercluster plane, and SGW chain denotes the superclusters in a chain in the foreground of the Sloan 
Great Wall. 



inner structure of superclusters in simulations from the present 
time to a distant future (from a = 1 to a = 100, a is the expan- 
sion factor). In their study superclusters were defined as high- 
mass bound objects, and the supercluster shape was approxi- 
mated by triaxial ellipses. To analyse the substructure of super- 
clusters they used the mult i plicity function of clusters in super- 
clusters. |^ra^a^Meloietal](|2009|) showed that superclusters are 
elongated, prolate structures, there are no thin pancakes among 
them. Future superclusters are typically much more spherical 
than present-day superclusters. Presently, the superclusters con- 
tain a large number of clusters, which may merge into a single 
cluster in the far future, i.e., multispiders and multibranching fil- 
aments may evolve into simple spiders and filaments. 



Comparisons of the properties of rich and poor superclusters 
have revealed several differences between them. The mean and 
maximum number densities of galaxies in rich superclusters are 
higher than in poor superclusters. Rich sup erclusters are mor e 
asymmetrical than poor superclusters (Einasto et al. 2007^ 
Rich superclusters contain high-density cores (Einasto et all 
l2007cl) . The fraction of rich clusters and X-ray cl usters in rich 
superclusters is larger than in poor superclusters (feinasto et alj 
2001; Luparello et al. 201 1), and the core regions of the richest 
superclusters may conta in merging X-ray clusters (Bardel li et al.l 
120001; iRose etani2002h . However, we still lack a detailed anal- 
ysis of whether the differences between the properties of the 
galaxy and group content of rich and poor superclusters are re- 
lated also to the differences in morphology of superclusters. 
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In lEinasto et alJ (l2008h we compared the properties of the 
two richest superclusters from the 2dF Galaxy Redshift Survey, 
the superclusters SCI 126 (SCI 061 in the present study), and 
the Sculptor supercluster (SCI 9 in EOl). We used Minkowski 
functionals to quantify the fine structure of these superclusters as 
traced by different galaxy populations. Our calculations showed 
that in the supercluster SCI 126 the population of red, early 
type galaxies is more clumpy than the population of blue, late 
type galaxies, especially in the outskirts of the supercluster. In 
contrast, in the supercluster SCI 9 the dumpiness of galaxies 
of different type is quite similar in its outskirts. In the core of 
the supercluster SCI 9 the dumpiness of blue, late type galax- 
ies is larger than the dumpiness of red, early type galaxies. 
In the supercluster SCI 1 1 1 in the SGW t he dumpiness of re d 
galaxies is larger than that of blue galaxies (lEinasto et al.ll20lTI) . 
Morphologically the supercluster SCI 126 resembles a multi- 
branching filament, while the Sculptor supercluster and SCI 111 
resemble a multispider. We need to study the morphology and 
galaxy populations of a larger sample of superclusters to find 
out whether the differences between galaxy populations in su- 
per clusters are also related to t heir different morphology. 

lAra^on-Calvo et aH (l2010h recently studied the structure and 
morphologies of elements that define the cosmic web with the 
Multiscale Morphology Filter and data from ACDM simula- 
tions. The authors found several typical morphologies for fila- 
ments, which they describe as line, grid, star, and complex fil- 
aments. We can compare that with our classification of super- 
cluster morphologies. Approximately, line filaments are com- 
parable with our simple filaments. Star filaments can be com- 
pared with simple spiders. Grid and complex filaments may cor- 
respond either to multibranching filaments or to multispiders, 
although complex filaments are more similar to multispiders and 
grid filaments to multibranching filaments. This shows that the 
morphology of observed and simulated superclusters is, in gen- 
eral, similar, as we found earlier from much smaller samples 
of observed and simulated superclusters. The biggest exception 
is the supercluster SCI 061, a very rich and high-density multi- 
branching filament. No other supercluster with such a morphol- 
ogy has been fo und yet either in simulations in observations (see 
also iGott et al . 2008). Another difference between the observed 
and simulated superclusters, as quantified by Minkowski func- 
tionals and shapefinders, is fine structure of superclusters, delin - 
eated by galaxies of different luminosity (lEinasto et al.ll2007dh . 
The dumpiness of observed superclusters for galaxies of differ- 
ent luminosity has a much larger scatter than that of simulated 
superclusters. Simulations do not yet explain all the features of 
observed superclusters. 

6. Summary and conclusions 

We have presented an analysis of the large-scale distribution 
and morphology of superclusters from the SDSS DR7. While 
the overall shape of superclusters has been analysed earlier in 
several studies, our paper is the first in which the inner mor- 
phology of a large sample of observed superclusters is studied in 
detail. We used multidimensional normal mixture modelling to 
divide superclusters into two sets according to their physical and 
morphological properties. We present 2D and 3D distributions of 
galaxies and rich groups in superclusters, the dumpiness curve 
(the fourth Minkowski functional Vi - mf relation), as well as 
the morphological signature K1-K2 for each supercluster in our 
sample . 

The superclusters were selected to contain at least 300 ob- 
served galaxies. As we showed, this does not mean that all su- 



perclusters in our sample are rich - about 1/4 of our superclusters 
contain one rich cluster or group only (three superclusters do not 
contain any rich group or cluster). 

Summarising, our study showed that: 

1) Most superclusters contain Abell clusters, about 1/3 of them 
are X-ray clusters. 

2) The large-scale distribution of the superclusters is very in- 
homogeneous. Rich superclusters in the sample form three 
chains, the Sloan Great Wall is the richest of them. 

3) Almost all superclusters under study are elongated and have 
filamentarities that are larger than than their planarities. 
More elongated superclusters are also more luminous, have 
larger diameters and contain a larger number of rich clusters. 
The values of the fourth Minkowski functional V3 show that 
they also have a more complicated inner morphology than 
less elongated superclusters. 

4) The morphological analysis shows a large variety of mor- 
phologies among superclusters. The fine structure of super- 
clusters can be described with four main types of morphol- 
ogy: spiders, multispiders, filaments, and multibranching fil- 
aments. Often a supercluster has an intermediate morphol- 
ogy between these main types. Superclusters with a sim- 
ilar shape parameter may have a different fine structure. 
Consequently neither the shape parameter or the number of 
rich clusters in a supercluster alone are sufficient to describe 
the morphology of superclusters. 

Our results on the morphology of superclusters can be used 
to compare the properties of local and high-redshift superclus- 
ters. Few superclusters at very high redshifts have alrea dy been 
discovered in deep, wide-field imaging surveys (N akata et alJ 
2005; Swinbank ettn i2007l: lOal et al.ll2008l:lTanaka et al. 200i; 
Schirmer et al .. iloTu). Deep surveys like the ALHAMBRA 
project (Moles et al I2OO8) will provide us with data about (pos- 
sible) very distant superclusters; we can analyse their structure 
and compare that with the local superclusters, using morpholog- 
ical methods. 

Our study does not give a definite answer to the question 
about the possible connection between the morphology of super- 
clusters and their large-scale distribution. Also, there are super- 
clusters in our sample that can be described as filaments or multi- 
branching filaments, but none of them is as rich and has such an 
overall high density as the supercluster SCI 061. Even the rich- 
est supercluster with a multispider morphology in our sample, 
the supercluster SCI 094, is not as ric h as the very rich Sculptor 
supercluster in lEinasto et al.l (l2007dl) . This shows the need for 
a larger sample of superclusters to understand the morphologi- 
cal variety of superclusters, and to study the possible connection 
between the large-scale distribution of superclusters and their 
morphology. 

Different morphologies of superclusters suggest that their 
evolution has been different. To understand the formation and 
evolution of superclusters of different morphology better, we 
need to study the properties and evolution of superclusters in 
simulations. The morphology of superclusters and its evolu- 
tion may be one of the factors to distinguish between different 
cosmological models (Kolokotronis et al. 2002; Hoffman et al] 
2007). Especially interesting is the supercluster SCI 061 in the 
Sloan Great Wall. Up to now simulations hav e not been able 
to model its morphology (lEinasto et alJ l2007dl and references 
therein). New simulations with larger volumes are needed to 
study the morphology of superclusters and the evolution of the 
morphology of simulated superclusters, to understand the rea- 
sons for the exceptional morphology of the supercluster SCI 06 1 . 
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In addition, very rich superclusters are rare; another reason to 
use simulations for large volumes comes from the demand to in- 
clude as large a variety of superclusters in the simulation volume 
as possible. 
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Appendix A: Luminosity density field and 
superclusters 

To calculate the luminosity density field, we must first calculate 
the luminosities of groups. In flux-limited samples galaxies out- 
side the observational window remain unobserved, and we have 
also to take into account the luminosities of these galaxies as 
well. For that, we multiply the observed galaxy luminosities by 
the luminosity weight W^. The distance-dependent weight factor 
Wd was calculated as follows: 



J^' Ln(L)dL' 



(A.l) 



where L12 = are the luminosity limits of the 

observational window at a distance d, corresponding to the ab- 
solute magnitude limits of the sur vey Mi and M2; we took 
Mq - 4.64 mag in the r-band (Blant on & Roweisll2007l) . Owing 
to their peculiar velocities, the distances of galaxies are some- 
what uncertain; if the galaxy belongs to a group, we used the 
group distance to determine the weight factor. 

To calculate a luminosity density field, we converted the 
spatial positions of galaxies r, and their luminosities L, into 
spatial (lurninosi ty) densities. For that we use kernel densities 
(ISilvermanll986l) : 



p(r) = ^K(r-rr.,a)Li, 



(A.2) 



where the sum is over all galaxies, and K (r; a) is a kernel func- 
tion of a width a. Good kernels for calculating densities on a 
spatial grid are generated by box splines B j. Box splines are lo- 
cal and they are interpolating on a grid; 



We defined the (one-dimensional) B3 box spline kernel K^^^ of 
the width a as 



K^il\x;a,6) = Bi{x/a){6/a), 



(A.5) 



where 6 is the grid step. This kernel differs from zero only in the 
interval x e [-2a, 2a]; it is close to a Gaussian with cr - I in 
the re gion x e [-a, a], so its effective width is 2a (see, e.g.. lSaai] 
l2009j) . The kernel exactly preserves the interpolation property 
for all values of a and 6, where the ratio a/6 is an integer. (This 
kernel can be used also if this ratio is not an integer, and a » 5; 
the kernel sums to 1 in this case, too, with a very small error). 
This means that if we apply this kernel to points on a one- 
dimensional grid, the sum of the densities over the grid is exactly 

The three-dimensional kernel /T^' is given by the direct 
product of three one-dimensional kernels: 



K''^\r; a, 6) = kI^\x; a, S)Kl^\y; a, 6)Kf\z; a, 6), 



(A.6) 



where r = {x,y,z]- Although this is a direct product, it is 
isotropic to a good degree (Saar 2009). 

The densities were calculated on a cartesian grid based on 
the SDSS 77, A coordinate system, because it allowed the most 
efficient fit of the galaxy sample cone into a brick. Using the 
rms velocity cr,,, translated into distance, and the rms projected 
radius cr,. from the group catalogue (TIO), we suppressed the 
cluster finger redshift distortions. We divided the radial distances 
between the group galaxies and the group centre by the ratio of 
the rms sizes of the group finger: 



%al,f 



^group ^" (^gal,i 



(A.7) 



This removes the smudging effect the fingers have on the density 
field. We used an 1 /i ' Mpc step grid and chose the kernel width 
a = 8 /i ' Mpc. This kernel differs from zero within the radius 
16 Mpc, but significantly so only inside the 8 /i"' Mpc radius. 

Before extracting superclusters w e applied the DR7 
mask constructed by P. Amalte-Mur (iMartmez et al.l l2009l 
iLiivamagi et al.ll2010h to the density field and converted densi- 
ties into units of mean density. The mean density is defined as 
the average over all pixel values inside the mask. The mask is 
designed to follow the edges of the survey and the galaxy distri- 
bution inside the mask is assumed to be homogeneous. 



Appendix B: Minkowski functionals and 
shapefinders 

For a given surface the four Minkowski functionals (from the 
first to the fourth) are proportional to the enclosed volume V, the 
area of the surface S , the integrated mean curvature C, and the 
integrated Gaussian curvature x- Consider an excursion set F^^ 
of a field 0(x) (the set of all points where the density is higher 
than a given limit, 0(x > (po)). Then, the first Minkowski func- 
tional (the volume functional) is the volume of this region (the 
excursion set): 



(A.3) VoiM 



d^x . 



(B.l) 



for any x and a small number of indices that give non-zero values 
for B j{x). We used the popular B3 spline function: 



The second Minkowski functional is proportional to the surface 
area of the boundary of the excursion set: 



B,ix) = ^ (1^ - 2\' - 4|x - 1|3 + 6\x\' - 41;. + Ip + I:. + 21^) .(A.4)^i(^o) = ^ d5(x) 



(B.2) 
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(but it is not the area itself, notice the constant). The third 
Minkowski functional is proportional to the integrated mean cur- 
vature C of the boundary: 



V2(cf>o) 



6^ Xf.„ 



1 

(X) 



dS(x) , 



(B.3) 



where R\(x) and Riix) are the principal radii of curvature of the 
boundary. The fourth Minkowski functional is proportional to 
the integrated Gaussian curvature (the Euler characteristic) of 
the boundary: 



Vi(cf>o) 



1 



47r JsF,„ Ri{x)R2ix) 



dS(x) 



(B.4) 



At high (low) densities this functional gives us the number of 
isolated clumps (void bubbles) in the sample jMartmez & Saarl 
[2002; Saaret al. 2007): 



^3 -^clumps -^bubbles -^tunnels ■ 



(B.5) 



As the argument labelling the isodensity surfaces, we chose 
the (excluded) mass fraction mf - the ratio of the mass in the 
regions with the density lower than the density at the surface, 
to the total mass of the supercluster When this ratio runs from 
to 1, the iso-surfaces move from the outer limiting boundary 
into the centre of the supercluster, i.e., the fraction mf - cor- 
responds to the whole supercluster, and mf = 1 - to its highest 
density peak. 

We used directly only the fourth Minkowski functional 
in this paper; t he other functionals were used to calculate 
the sh apefinders (ISahni et al.lll998l: IShandarin et"aD '2004'; 'Saar' 
12009 ). The shapefinders are defined as a set of combinations 
of Minkowski functionals: Hi - 3V/S (thickness), H2 - S/C 
(width), and Hj - C/4n (length). The shapefinders have dimen- 
sions of length and are normalized to give //, = R for a sphere of 
radius R. For smooth (ellipsoidal) surfaces, the shapefinders //, 
follow the inequalities Hi < H2 < Hj,. Oblate ellipsoids (pan- 
cakes) are characterised by Hi « H2 ~ H^, while prolate ellip- 
soids (filaments) are described hy Hi ^ H2 « H^. 

[Sahni et al. (1998) also defined two dimension- 
less shapefinders Ki (planarity) and K2 (filamentarity): 
Ki = iH2 - Hi)l{H2 + Hi) and K2 = (H^ - H2)/{Hi + H2). 

In the (/Ti, /r2)-plane filaments are located near the /Ta-axis, 
pancakes near the Ki -axis, and ribbons along the diameter, con- 
necting the spheres at the origin with the ideal ribbon at (1 , 1). In 
lEinasto et al.l (l2007dh we calculated typical morphological sig- 
natures of a series of empirical models that serve as morpho- 
logical templates to compare with the characteristic curves for 
superclusters in the {Ki , /r2)-plane. 



Appendix C: 3D figures of superclusters 

Here we show the 3D distribution of groups in our superclusters. 
Figure ICTI is 3D presentation of Fig.|3]and shows the distribu- 
tion of groups with at least twelve member galaxies in our su- 
perclusters in the cartesian coordinates defined in Section l4n in 
units of /; ' Mpc. The filled circles denote groups with at least 
50 member galaxies, empty circles denote groups with 30^9 
member galaxies and crosses denote groups with 12-29 mem- 
ber galaxies. The numbers are ID numbers of superclusters with 
at least 950 member galaxies. Fig. lC.2l IC.4I and IC.4l show 3D 
distribution of groups in superclusters. In these figures we plot 
rich groups with at least 30 member galaxies with filled circles. 



groups with 2-29 member galaxies with empty circles. We do 
not show the figures for superclusters, for which the value of the 
fourth Minkowski functional Vj, - \ over the whole mass frac- 
tion interval. 
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Fig. C.l. Distribution of groups with at least 12 member galaxies 
in our superclusters in the cartesian coordinates as in Fig. [3] in 
units of Mpc. Filled circles denote the groups with at least 
50 member galaxies, empty circles denote the groups with 30- 
49 member galaxies and crosses denote the groups with 12-29 
member galaxies. The numbers show the ID's (Table [1] column 
1) of very rich superclusters from Sect. 14.31 



Appendix D: The distribution of galaxies and rich 
groups in the sky and the fourth Minkowski 
functional and morphological signature for 
superclusters with less than 950 member 
galaxies 

In this section we present the distribution of galaxies and rich 
groups/clusters in the sky as well as the fourth Minkowski func- 
tionals and morphological signatures for superclusters with less 
than 950 member galaxies (Fig.|ET] |D21 and ID3] ). In all fig- 
ures the left panels show the distribution of galaxies in the sky 
(dots). Circles in these panels mark the location of groups with 
at least 30 member galaxies. The size of a circle is proportional 
to the size of a group in the sky. The middle panels show the 
fourth Minkowski functional vs. the mass fraction mf, and 
the right panels show the shapefinders Ki (planarity) and K2 (fil- 
amentarity) for the supercluster (the morphological signature). 
Filled circles in the right panels mark the value of the mass frac- 
tion mf - 0.7. In the right panels the mass fraction increases 
anti-clockwise along the curves. As in the previous section, we 
do not show the figures for superclusters, for which the value of 
the fourth Minkowski functional Yt, - \ over the whole mass 
fraction interval. 
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Fig. C.2. 3D figures of the very rich superclusters described in 
Sect. 14.31 Filled circles show the location of rich groups and 
clusters with at least 30 member galaxies, empty circles show 
poorer groups. We plot the right ascension (in degrees), decli- 
nation (in degrees), and distance (in Mpc//i) of groups. Plots of 
superclusters are given in the same order as they are presented 
in the text. From top to bottom: left; the supercluster SCI 061, 
right: SCI 094, left; the supercluster SCI 001, right; SCI Oil, 
left; the supercluster SCI 024, right; SCI 055, left; the superclus- 
ter SCI 336, right; SCI 350. 































58 157 156 


55 154 1 





fl77.0 176.5 176.0 175.S 175.0 174.5 174.0 

















Co ^ 




*°°o 

















-^12 211 210 209 208 207 206 205 
R.A. 




2 210 20S 206 



^50 258 257 256 255 254 253 252 
R.A. 



Fig. C.3. 3D figures of superclusters with less than 950 mem- 
ber galaxies. Panels and notations in this figure are the same as 
in Fig. IC.2I From top to bottom; left; the supercluster SCI 038, 
right; SCI 064, left; the supercluster SCI 087, right; SCI 136, 
left; the supercluster SCI 152, right; SCI 189, left; the super- 
cluster SCI 198, right; SCI 223, left; the supercluster SCI 228, 
right; SCI 317, left; the supercluster SCI 332, right; SCI 349, 
left; the supercluster SCI 351, right; SCI 362, left; the superclus- 
ter SCI 366, right; SCI 376. 
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Fig. C.4. 3D figures of superclusters with less than 950 member 
galaxies (continued). Panels and notations in this figure are the 
same as in Fig. IC.2I From top to bottom: left: the supercluster 
SCI 474, right: SCI 512, left: the supercluster SCI 525, right: 
SCI 530, left: the supercluster SCI 779, right: SCI 827. 
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Fig.D.l. Left panels show the distribution of galaxies (dots) in 
the sky. Circles mark the location of groups with at least 30 
member galaxies, and the size of a circle is proportional to the 
size of a group in the sky. The middle panels show the fourth 
Minkowski functional vs. the mass fraction mf, and the right 
panels show the shapefinders K\ (planarity) and K2 (filamer^ 
tarity) (the morphological signature). In the right panels filled 
circles mark the value of the mass fraction mf - 0.7; the mass 
fraction increases anti-clockwise along the curves. From top to 
bottom: the superclusters SCI 038, SCI 064, SCI 087, SCI 136, 
SCI 152, SCI 189, SCI 198, and SCI 223. 
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Fig. D.3. Panels in this figure are the same as in Fig. ID. II From 
top to bottom: the superclusters SCI 474, SCI 512, SCI 525, 
SCI 530, SCI 779, and SCI 827. 



Fig.D.2. Panels in this figure are the same as in Fig. ID. II From 
top to bottom: the superclusters SCI 228, SCI 317, SCI 332, 
SCI 349, SCI 351, SCI 362, SCI 366, and SCI 376. 
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